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Abstract 

Kinetic equations are analyzed for thermal degradation of polymers. The governing 
relations are based on the fragmentation-annihilation concept. Explicit solutions to these 
equations are derived in two particular cases of interest. For arbitrary values of adjustable 
parameters, the evolution of the number-average and mass-average molecular weights of 
polymers is analyzed numerically. Good agreement is demonstrated between the results of 
numerical simulation and experimental data. It is revealed that the model can correctly 
predict observations in thermo-gravimetric tests when its parameters are determined by 
matching experimental data for the decrease in molecular weight with exposure time. 
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1 Introduction 



This paper is concerned with the kinetics of thermal degradation of polymers. This subject has 
attracted substantial attention in the past half a century both among the specialist in theoretical 
physics and chemical engineering |2|. This may be explained by two reasons: (i) scission 
(fragmentation) of macromolecules driven by thermal fluctuations at elevated temperatures 
provides a good example for the analysis of population dynamics in complex systems, and (ii) 
analysis of the degradation process becomes more and more important due to an increase in the 
range of temperatures for engineering applications, recycling of post-consumer plastic waste, 
as well as the use of polymers as biological implants and matrices for drug delivery, where 
depolymerization is an inevitable process affecting the life-time of an article. 

Although the fact that binary scission of chains is the main mechanism for thermal degra- 
dation of polymers is widely accepted, the applicability of the fragmentation concept to the 
analysis of experimental data is rather limited. This may be attributed to the fact that ki- 
netic equations for the "linear" fragmentation process involve only one adjustable parameter 
jHj, the number too small to provide a reasonable approximation for observations. The latter 
implies that generalizations of conventional equations for the kinetics of fragmentation become 
a meaningful subject for investigation. 

The following ways for "refinement" of the fragmentation concept may be mentioned: 

1. An account for the influence of chains' length and inhomogeneity of chains (the fact that 
the probability of a breakage event depends on the position of a bond along the backbone 
of a chain) on the rate of scission |U El E] • 

2. An introduction of two kinds of macromolecules in a network (with strong and weak 
bonds) that have different rates of fragmentation [TT H ITT | ITS]. 

3. An increase in the number of species in the population balance laws by accounting for 
interactions between chains and free radicals and the effect of the concentration of free 
radicals on the scission rate [THl H3J [Tj3 HBJ H3 HHj . 

Although these approaches seem rather attractive from the theoretical viewpoint, they are 
grounded on physical assumptions that are hard to be verified in experiments, on the one hand, 
and they result in overcomplicated explicit solutions (when the latter can be developed), on 
the other. 

Another way to make the fragmentation concept more flexible for matching observations is to 
take into account some physical processes at the micro-level that accompany thermally-induced 
scission of macromolecules. Two candidates for this role appear to be natural: (i) aggregation 
(recombination) of broken chains [TTHl2THl2~Tl l22|. an d (ii) annihilation of fragmentation products 
(either in the form of creation of inert species that do not take part in further fragmentation 
[2ID I2U 121] or as diffusion of small-size fragmentation products and their subsequent evaporation 
through the surface of a specimen 123 I2S ) • 

In this study, we adopt the latter approach and accept the conventional hypothesis that the 
diffusivity of detached end- and side-groups at the degradation temperature is so large that 
the kinetics of diffusion may be disregarded. This allows the number of adjustable parameters 
in the governing equations to be reduced to four, the number that gives an opportunity to 
analyze the influence of all these quantities on the evolution of molecular weight of polymers. 
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Analysis of the effect of diffusivity of oligomers on the degradation kinetics will be the subject 
of a subsequent paper. 

Unlike previous studies, which focused mainly on the scaling solutions to the fragmentation- 
annihilation equations (the latter provide valuable information about the properties of the 
degradation process at small and large times), we concentrate on the kinetics of polymer degra- 
dation within the time-scale of conventional tests. Another difference between the present work 
and previous publications is that we deal with changes in the number-average and mass-average 
molecular weights, the standard quantities measured in experiments. This allows our analytical 
and numerical results to be compared with available observations. 

The objective of this study is two-fold: 

• To report some analytical and numerical solutions to the fragmentation-annihilation equa- 
tions. 

• To establish restrictions on the material constants in the governing relations imposed by 
the time-temperature and time-molecular weight superposition principles, and to find 
adjustable parameters in the kinetic equations by matching observations. 

We believe that the knowledge of these quantities (or, at least, their orders of magnitude 
and mutual relations between different constants) may be helpful for further analytical and 
numerical investigation of the degradation process. 

In fitting experimental data, we concentrate on thermal degradation of polystyrene (PS), 
poly(a-methylstyrene) (PAMS) and poly(L-lactide) (PLLA). The choice of these polymers may 
be explained by their wide use in industrial applications, including food service packaging (PS 
and PAMS) and scaffolds for transplanted organs (PLLA), the areas, where the degradation 
process is of the highest importance. Another (merely technical) reason for this choice is that 
observations are available on these polymers with low polydispersity indices (less than two), 
which implies that no additional parameters associated with the initial distribution of chains 
should be introduced into the model. 

The paper is organized as follows. We begin with the formulation of kinetic equations in 
Section 2. Some explicit solutions to these equations are derived in Section 3. The effect of 
material constants on the evolution of molecular weight of polymers subjected to thermal degra- 
dation is studied numerically in Section 4. The time-molecular weight and time-temperature 
superposition principles are discussed in Sections 5 and 6, respectively. In Section 7, we demon- 
strate that the model can adequately predict observations in thermo-gravimetric tests. Some 
concluding remarks are formulated in Section 8. 

2 Kinetic equations 

In this section, kinetic equations are formulated for the analysis of thermal degradation of 
polymers based on the fragmentation-annihilation concept. 

2.1 Binary scission of chains 

Denote by N(t) the number of macro molecules per unit mass of a polymer network at an 
arbitrary instant t > 0. Following common practice, we treat chains as sequences of segments 
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connected by bonds. Each chain in the network is entirely characterized by the numbers of 
segments k (k — 1, 2, . . .). Denote by N k (t) is the number of chains (per unit mass) at time t 
containing k segments. The functions N k (t) obey the conservation law 



N(t) = J2N k (t). (1) 



k=l 



Binary scission (fragmentation) of chains is described by the reactions 

N k - Nt + (l = l,...,k-l). 

Denote by 7 the rate of scission (the number of scission events per bond between segments per 
unit time). Assuming 7 = 70 to be independent of the number of segments k, we arrive at the 
kinetic equations for the functions N k (t) 



00 

(t) = - 7o (fc - l)N k (t) + 2 7o £ Nj(t). (2) 



dNk 

dt j=k+ i 

The coefficient k — 1 in the first term describes the number of possible scission events in a chain 
containing k segments. The coefficient "2" before the sum in Eq. (2) indicates that there are 
two opportunities ("left" and "right") to obtain a chain with k segments after scission of a 
macromolecule with a larger number of segments. 

Equation (2) can be easily generalized by assuming the rate of scission 7 to be a function 
of chains' length. Conventionally, the power-law dependence is adopted to account for this 
dependence 

7 (fc) = 7o* fl , (3) 

where 70 and a are material parameters. The kinetic equations for fragmentation of chains with 
the length-dependent rate of scission (3) read 

°-±{t) = - lo k a (k - l)N k (t) + 2 7o £ 3 a Nj(t). (4) 

j=k+l 



2.2 Annihilation of chains 

Thermal fluctuations in a network induce not only binary scission of macromolecules, but also 
detachment of end- and side-groups from polymer chains. As these groups are rather small, 
they have relatively large diffusivity and can easily leave a specimen. A decrease in a sample's 
mass with time driven by separation of end- and side-groups and their desorption is treated as 
their annihilation. 

To make the model tractable from the mathematical standpoint, we suppose that detach- 
ment of small groups within the interval [t, t + dt] may be thought of as transformation of a 
chain with k segments into a chain with k — 1 segment. Denote by T the ratio of the number 
of chains that lose a segment per unit time to the entire number of chains. Assuming the an- 
nihilation rate to grow with the number of segments in a chain following the power law similar 
to Eq. (3), 

r(k) = r k b , (5) 



4 



where T and b are material parameters, we arrive at the kinetic equations 

dN ki 



dt 



-(t) = - lo k a (k - l)N k (t) + 2 7o + r o (k + l) b N k+1 (t) - k b N k (t) . (6) 



j=k+i 



The number-average molecular weight M n and the mass-average molecular weight M w of a 
polymer read 



(7) 



The objective of this work is to analyze changes in M n (t) and M w (t) with time t, when the 
functions N k (t) are governed by Eq. (6). 



2.3 Transformation of the kinetic equations 

Introducing the concentrations of chains with k segments, 

*M = (8) 

where iVo = N(0) is the total number of chains at the initial instant t = 0, we present Eq. (6) 
in the form 

dni 00 r n 

-±(t) = - lQ k a (k-l)n k (t) + 2 l0 j a nj(t)+r [{k + l) b n k+1 {t)-k b n k (t)]. (9) 

di j=k+l 

Following common practice, it is convenient to suppose that the number of segments in a chain 
is large compared to unity and to replace the discrete index k in Eq. (9) by a continuous 
argument x. This results in the integro-differential equation for the function n(t, x) 

— {t,x) = - lo x a+1 n{t,x)+2 1 oJ x y a n{t,y)dy + r — (x b n{t,x)). (10) 

The initial condition for Eq. (10) reads 

n(0,x) = n (x), (11) 

where no(x) is a given function. We do not formulate boundary conditions for the function 
n(t, x), but assume that this function does not grow very strongly at x = and decays rapidly 
at x — > oo in the sense that the integrals M m exist, 

POO 

M m {t)= x m n{t,x)dx (m = 0,1,2,...). (12) 



3 Explicit solutions 

At T = 0, i.e. when the annihilation process is suppressed, an analytical solution of Eq. (10) 
was found in [5] for an arbitrary integer a > 0. However, an appropriate solution for a > is 
so complicated that it can say practically nothing about the behavior of the function n(t, x) 
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and its moments M m (t). This implies that we confine ourselves to the case a = in the search 
for explicit solutions of Eq. (10). Following [5], we assume the parameter b in Eq. (10) to be 
an integer as well, b — I (I — 1, 2, . . .), and re-write Eq. (10) in the form 

OtI /*oo 0* / \ 

-r^(t,x) = -j Q n{t,x) + 270 y n(t,y)dy + T — (xn(t,x)). (13) 
Integrating Eq. (13) over x m (m — 0, 1, . . .) and using notation (12), we find that 

^(£) = -7oM m+1 (t) + 2 7o y o x m dxj x n(t,y)dy + T j o x m — (x l n{t,x))dx. (14) 
The first integral in Eq. (14) is transformed by changing the order of integration, 

■ poo poo py | poo | 

x m dx / n{t,y)dy = / n(t,y)dy / x m dx = — / y m+1 n(t,y)dy = — M m+1 (t). 

Jx Jo Jo m + 1 jo m + 1 

The other integral is calculated by integration by parts, 

/ x m — (xn(t,x))dx = -m x rM n(t,x)dx = -mM m+w (i). 
jo ox v ' Jo 

Substitution of these expressions into Eq. (14) results in the differential equation 

-^(t) = -7o (l 2 —)M m+l (t) - r mM m+H (f). (15) 

at V m + 1 / 

Our purpose now is to find explicit solutions of Eq. (15) for / = 1 and / = 2. According to Eq. 
(5), the former case corresponds to scission and annihilation of end-groups exclusively, whereas 
the latter case describes annihilation of side-groups homogeneously distributed along a chain's 
backbone. 



3.1 The case / = 1 

Setting I — 1 in Eq. (15), we obtain 

dM m 
dt 

Introducing the notation 



(t) = -T mM m (t) - 7o (l —]M m+1 (t). (16) 

V 171+1/ 



T m (r) = M m (t)exp(T mt), t=( exp(-T Q s)ds, (17) 

Jo 



we transform Eq. (16) as follows: 

dT m/ , / 2 



dr 



(T) = -7o(l-^ n -V m+1 (r). (18) 
V m + 1/ 



Comparison of Eqs. (15) and (18) implies that in the special case / = 1, the kinetics of the 
fragmentation-annihilation process coincides with the kinetics of a mere fragmentation process 
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in the new time r for the new moments T m . In particular, putting m = and m — 1 in Eq. 
(18), we find that 



dT 
dr 



7o^ 



i- 



dTi 
dr 



0. 



The solutions of these equations are given by 

T (r) = T (0) + 7 oTi(0)r, T x {r) = T\(0). 

The moment T 2 (r) is determined by the formula 

2 



r 2 (r) 



(7ot) 



n (7or)-T (0) + 71(0)7,^ 



(19) 



(20) 



where no(^) is the Laplace transform of the function Uq(x). A detailed derivation of Eq. (20) 
is given in Appendix. It follows from Eqs. (17), (19) and (20) that 



M (t) 
M 1 {t) 
M 2 (t) 



M (0) + M 1 (0)^(l-exp(-r t) 
t o v 



M!(0)exp(-r t 



TV exp(-2r t) 



7o 1 (l-exp(-r t)) 2 



^o(^(l 

VJ - 



exp(-iy 



-M (0) + Mx(0)^(l - exp(-r t 
J- o v 



(21) 



Substituting Eq. (21) into Eq. (7), where the sums are replaced by integrals, we arrive at 
explicit expressions for the number-average and mass-average molecular weights as functions of 
time t. 



3.2 The case / = 2 

Setting I = 2 in Eq. (15), we obtain 

^W = -7o(^ + ^mW m+lW . (22) 
at V m + 1 70 / 

Assuming the initial network of chains to be monodisperse, 

n {x) = 5{x-L), (23) 

where L is the initial length of chains, we solve Eq. (22) by Charlesby's method [HO]- As 
Eq. (22) is linear with respect to the unknown function n(t,x), appropriate formulas for the 
moments M m (t) corresponding to an arbitrary initial condition no(x) are developed by the 
superposition method. 

The mth moment M m (t) is expanded into the Taylor series in time, 

^ (t) = f^, ,24) 

k=0 
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where M^(0) stands for the kth derivative at the point t — 0. Substitution of Eq. (24) into 
Eq. (22) implies that for an arbitrary k > 1, 



fc-i 



^(0) = (- 70 )* J] (^4^ + ^(m + j))M m+k (0). 
It follows from Eqs. (12) and (23) that 

M m (0) = L m . 

Substitution of expressions (25) and (26) into Eq. (24) results in 



M m (t) 



l + EAnfc(-70^)' 

k=l 



(25) 

(26) 
(27) 



where 



Imfc 



1 ^Vm + j - 1 



;,n 



+ — (m + j) 
m + j + 1 7o 



Introducing the new variable f = j + 1 and omitting the prime, we obtain 



1 fm + j - 2 T . 

— — + — (m + j-1 . 

m + j 7o 



(28) 



It follows from Eq. (27) that for an arbitrary initial condition no(x), the moments M m (t) are 
given by 



roc r °° 

M m (t) = / n {x)x m 1 + E A mk {- lo tx) k 

CO 

= M m (0) + J2 A rnkM m+k (0)(-^t) k - 



dx 



(29) 



k=i 



Although Eq. (29) provides explicit expressions for the moments M m (t), it is not rather con- 
venient for the numerical analysis, because the series converges slowly. Formulas (28) and (29) 
are helpful, however, for the evaluation of changes in M m (t) at small times, when 7 t 1. 
Neglecting terms beyond the first order of smallness in Eq. (29), we find that 

M (t) = M (0) + M 1 (0) 7o i, M 1 (t) = Mx(0) - M 2 (0)r t, 



M 2 (t) = M 2 (0) - M 3 (0)(-^ + 2r )*. 



(30) 



According to Eq. (30), changes in the moments M (t) and Mi(t) are governed by two different 
processes: an increase in M (t) is driven by fragmentation of chains, whereas a decrease in 
Mi(t) is induced by annihilation of end- and side-groups. Introducing the notation 



d(t) 



Mi(t) 



M n (t) 



d w {t) 



M w (t) 
M w (0) 



Mi(0)' w M n (0)' 

and using Eq. (7), we present Eq. (30) in the form 

d(t) = i - M w (o)iy, d Q (t) = i - (M n (o) 7o + M w (o)r )t, 



(31) 



dw{t) 



M w (o)r + M z (o)(^ + 2r ( 



t, 



(32) 



where the conventional notation is used M z (t) = M 3 (t)/M 2 (t). 
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4 Numerical analysis 



As explicit formulas for the moments M m (t) can be derived only for a limited set of values of 
the exponents a and b, it is of interest to analyze the evolution of M m with time t numerically 
for arbitrary values of these parameters. For this purpose, we integrate Eq. (9) numerically 
with the initial condition 

nk{0) = 8 k K, (33) 
where 6y denotes the Kronecker delta. Equation (33) corresponds to the monodisperse distri- 
bution of chains that contain K segments at the initial instant t = 0. To reduce the number of 
material parameters, we introduce the dimensionless time t' = 7ot and set 

V = -• (34) 
7o 

In the new notation, Eq. (9) reads (the prime is omitted for simplicity) 

d ^(t) = -k a {k-l)n k {t) + 2 £ j a nj{t)+ri[(k + l) b n k+1 (t)-k b n k (t)] (k = 1, 2, ... , K). 

at j=k+l 

(35) 

The limitation on the number of equations in Eq. (35) follows from the fact that if the maximal 
number of segments in a chain equals K at t = 0, no chains with higher number of segments 
can appear at t > due to the fragmentation-annihilation process. Equations (33) and (35) 
involve four material constants: a, b, r\ and K. In the numerical simulation, we fix K = 100 and 
study the effect of other parameters on the number- average and mass-average molecular weights 
determined by Eqs. (7) and (8). Integration of Eq. (35) is performed by the Runge-Kutta 
method with the step At = 1.0 • 10~ 6 . 

First, we set a = and 6 = 1 and study the effect of rj on the ratios d n (t) and d w (t) defined 
by Eq. (31). Figures 1 and 2 demonstrate that annihilation of end- and side-groups affects the 
degradation process noticeably when the rate of annihilation r exceeds the rate of scission 70 
by an order of magnitude. 

Afterwards, we fix 77 = 10.0, b = 1.0 and analyze the influence of the exponent a on the 
dimensionless ratios d n and d w . Figures 3 and 4 reveal that an increase in a results in a very 
strong decrease in d n and d w for any t > 0. When a « 1, the curves d n (t) and d w (t) steeply drop 
at small times and remain practically constant at larger times. These results make questionable 
the applicability of the analytical solution developed in (this solution corresponds to a = 2) 
for the description of the degradation process. 

Finally, we fix a = 0.2, 77 = 10.0 and assess the effect of the exponent b on d n (t) and d w (t). 
Figures 5 and 6 show that the influence of b is negligible when b G [0, 1). However, for b > 1, the 
effect of this parameter becomes substantial, and an increase in b causes a pronounced decrease 
in number- average and mass-average molecular weights. 

To ensure the accuracy of numerical simulation, we use three tests. First, we verify that at 
rj = 0, the first moment Mi remains independent of time [this conclusion follows from Eq. (15), 
where we set r = 0] . Secondly, we increase K by twice, decrease the rate of fragmentation 70 
by twice and check that the moments M n (t) and M w (t) remain unchanged. The latter implies 
that the results of numerical analysis are independent of our choice of K = 100. Finally, 
we perform simulation with a = and 6=1 and confirm that the numerical results for the 
moments M m (t) (m = 0, 1, 2) coincide with analytical solution (21). 
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5 Time— molecular weight superposition principle 



Experimental data for monodisperse polymers reveal that observations for the ratios of number- 
average and mass-average molecular weights, d n {t) and d w (t), obtained at various initial molec- 
ular weights and plotted in semi-logarithmic coordinates (versus the logarithm of time) may be 
superposed (with an acceptable level of accuracy) by shifts along the time axis. The construc- 
tion of a master-curve by shift of creep and relaxation curves measured at various temperatures 
is a conventional procedure in linear viscoelasticity of polymers [HI]. The validity of this op- 
eration is based on the time-temperature superposition principle, which asserts that only the 
characteristic creep and relaxation times are affected by temperature, while other parameters 
(elastic moduli and shapes of the relaxation spectra) remain temperature-independent. Our 
aim now is to establish an analogous principle for the effect of initial molecular weight of 
monodisperse polymers on the fragmentation-annihilation process. For this purpose, we return 
to integro-differential equation (10) and introduce the dimensionless variables 

t* = q, = — , (36) 

where G is the characteristic time for degradation. In the new notation, Eq. (10) reads 





-x a + 1 n(t*,x*) + 2 j(_ y y(t,,y,)dy, + r]L b - a - 2 — [xln(U,x,)) 



■ (37) 



It follows Eq. (37) that any solution n(t*,x*) is independent of the initial molecular weight L 
[this statement implies the time-molecular weight superposition principle] provided that 

b = a + 2. (38) 

Under condition (38), the shift factor A = 0/@ rcf , where B ref is the characteristic time corre- 
sponding to a reference molecular weight L ref , is determined by the formula 

logA=(a + l)log-^ (39) 

with log = log 10 . According to the time-molecular weight superposition principle, the number 
of adjustable parameters in the model may be substantially reduced: instead of four material 
constants, 70, T , a and b, we have only two parameters to be found, 70 and 77: the exponent a 
is uniquely determined from Eq. (39), whereas the exponent b is given by Eq. (38). 



5.1 Fitting of observations on polystyrene 

To find the values of 70 and 77 and to assess the influence of annihilation of end- and side-groups 
on the degradation process, we focus on observations reported by Madras et al. [11] on thermal 
degradation of monodisperse PS in a mineral oil at the temperature T = 275 °C. For a detailed 
description of specimens and the experimental procedure, the reader is referred to [IT] . The 
ratio d w of mass-average molecular weights is depicted versus exposure time t in Figure 7. The 
experimental data for M w (0) = 26 kg/mol are presented without changes. Observations for 
the other initial molecular weights (M w (0) = 12, 110, 210, 330 and 930 kg/mol) are shifted 
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along the time-axis by appropriate amounts A that are determined from the condition that all 
available data produce a smooth master-curve. 

The parameter A is plotted versus the ratio of initial molecular weights M w (0)/M^ ef (0) with 
M^ ef (0) = 26 kg/mol in Figure 8. The experimental data are approximated by the function 

^A = A + A 1 lo g ^^ y (40) 

where the coefficients A^ (k = 0, 1) are determined by the least-squares method. Figure 8 
demonstrates that Eq. (40) provides good matching of the observations with A ph and 
Ax pa 1. With reference to this result, we conclude from Eqs. (38), (39) and (40) that one 
should set a = and b = 2 in Eq. (35) in order to reproduce the experimental data plotted in 
Figure 7. 

The fact that « 1, which implies that the fragmentation rate 70 is proportional to the 
initial mass-average molecular weight M w (0), 

7o_ M w (0) 

7 -f M ref (0) > I ) 

has been revealed about 40 years ago [35] . It was also shown that Eq. (41) is valid when 
the initial molecular weight M w (0) is not very large: for polymers with ultra-high molecular 
weights, the fragmentation rate 70 grows like the square- root of the molecular weight [3*3*j . 

To find 70 and 77, we fix some intervals [0, 7 max ] and [0, r/ max ], where the "best-fit" parameters 
70 and r\ are assumed to be located, and divide these intervals into J subintervals by the points 
7W = 2A7, and r/ j) = jAi] (i,j — l,...,J — 1) with A7 = 7 max / J and Ai] = i] mSL J J. For any 
pair {7^,77^}, Eq. (35) with initial condition (33) is integrated numerically by the Runge- 
Kutta method with K = 100 and the time-step At = 0.1. This, relatively large, step is chosen 
because the fragmentation rates under consideration 70 are quite small (of order of 10 -5 ). The 
best-fit parameters 70 and r\ are determined from the condition of minimum of the function 

^ = E[C xp (U-C um (U] 2 , 

where the sum is calculated over all times t m at which observations are presented in Figure 7, 
d^ xp is the ratio of mass-average molecular weights measured in the tests, and d^ nm is given 
by Eqs. (7), (8) and (31). Figure 7 demonstrates fair agreement between the observations at 
various initial molecular weights and the results of numerical simulation with 70 = 1.4- 10 -5 and 
r\ = 1.0. Two conclusions may be drawn from the numerical analysis: (i) the model correctly 
describes the experimental data, and (ii) the best-fit value rj — 1.0 found by matching the 
observations on monodisperse polystyrene results in a rather weak influence of the annihilation 
process on the kinetics of thermal degradation (see Figures 1 and 2). Our purpose now is to 
demonstrate that this is not the case for other polymers. To show that annihilation of end- 
and side-groups noticeably affects the degradation process, we analyze two sets of observations 
on PAMS and PLLA obtained at various temperatures T. 



6 Time— temperature superposition principle 

We begin with the discussion of the effect of temperature T on the kinetics of thermal degrada- 
tion. To establish conditions on the material constants which ensure that the time-temperature 
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superposition principle holds for the ratios d n and d w , i.e. that the functions d n (t) and d w (t) 
measured at various temperatures and plotted in semi-logarithmic coordinates may be shifted 
along the time-axis to construct a master-curve ^j, we return to the dimensionless integro- 
differential equation (37). 

According to the Eyring theory the effect of temperature on the rates of thermally 
activated processes may be accounted for by the formulas 



7o = Y ei exp 



_2 
R 



Tf exp 



E-p / 1 1 1 



(42) 



where T is the absolute temperature, R is the universal gas constant, E n and are appropriate 
activation energies, and 7g ef and r^j are rates of fragmentation and annihilation at the reference 
temperature T ref . Equation (37) implies that the time-temperature superposition principle is 
valid, provided that the exponents a and b are independent of temperature, while the rates of 
fragmentation and annihilation, 70 and T , are affected by T in a similar fashion, 



a = a 



ref 



b = b 



ref 



En, = Et 



E. 



(43) 



The first two conditions in Eq. (43) seem quite natural, whereas the last equality imposes 
rather strong restrictions on the degradation process. In what follows, it will be shown that 
the latter condition is fulfilled for some polymers and is violated for others. 



6.1 Fitting of observations on poly(a— methylstyrene) 

We begin with the analysis of observations for thermo-oxidative degradation on PAMS with 
the initial mass-average molecular weight M w (0) = 9.0 kg/mol and the polydispersity index 
D = 2 reported by Sterling et al. These data demonstrate that the time-temperature 

superposition principle is satisfied for the degradation process with a high level of accuracy. 

In each test, a specimen was dissolved in 1,2,4-trichlorobenzene, heated to a fixed reaction 
temperature T in the range from 110 to 150 °C, and di-tert-butyl peroxide was added to a 
reaction vessel with a constant flow rate under stirring. Specimens for GPC (gel permeation 
chromatography) analysis were taken every 45 min. For a detailed description of the experi- 
mental procedure, see [To] . 

Changes in the ratio of number- average molecular weights d n with time t are depicted in 
Figure 9. The experimental data at T ref = 110 °C are presented without changes. Observations 
at the other temperatures (T = 120, 130, 140 and 150 °C) are shifted along the time-axis by 
appropriate amounts A that are determined from the condition that the experimental data 
produce a smooth master-curve. 

The parameter A is plotted in Figure 10 versus temperature T. The experimental data are 
approximated by the dependence [that follows from Eq. (42)] 

\nA = A -^ (44) 

with 

Tn E E 

A=-^j, A = 7, At = -. (45) 
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The coefficients (k = 0, 1) in Eq. (44) are determined by the least-squares method. Figure 
10 demonstrates that Eq. (44) ensures quite acceptable fit of the observations. According to 
Eq. (45) and Figure 10, the activation energy reads E = 86 kJ/mol. This value is in reasonable 
agreement with the activation energy E = 112 kJ/mol found by Brown and Wall [33] by 
applying a graphic method (whose accuracy is rather low due to the numerical estimation of 
derivatives) and in excellent accord with the activation energies recently determined for other 
polymers by using more sophisticated techniques (E = 88 kJ/mol for low-density polyethylene 
[T2*] and E = 98 kJ/mol for polypropylene JDj)- 

To reduce the number of adjustable parameters in Eq. (35), we set a = 0, which means that 
we presume the rate of fragmentation to be independent of chains' length. To find the quantities 
7o, b and 77, we use an algorithm similar to that employed in Section 5. We fix some intervals 
[0,7max], [0, 6 m ax] and [0, T] m3X \, where the "best-fit" parameters 70, b and r\ are assumed to be 
located, and divide these intervals into J subintervals by the points 7W = zA7, b^> = jAb and 



n 



kAr] k = 1, . . . , J - 1) with A7 = 7 max / J, Ab = 6 max / J and Ar] = r? max / J. For 
any triple {7^,6^,77^}, Eq. (35) with initial condition (33) is integrated numerically by the 
Runge-Kutta method with K = 100 and the time-step At = 0.1. The best-fit parameters 70, 
b and rj are determined from the condition of minimum of the function 



# = EK exp (U-< num (^ 



(46) 



where the sum is calculated over all times t m at which observations are presented in Figure 9, 
d£ xp is the ratio of number- average molecular weights measured in the tests, and d^ nm is given 
by Eqs. (7), (8) and (31). Figure 9 demonstrates good agreement between the observations at 
various temperatures T and the results of numerical simulation with 7 = 9.0 • 10~ 7 , b = 1.5 
and rj = 50.0. 

As the value of r] is rather large compared to unity, we draw a conclusion from Figures 
1, 2 and 9 that the effect of annihilation of end- and side-groups on the degradation kinetics 
of PAMS is substantial. The observation that the rate of annihilation for PAMS (rj = 50) 
exceeds that for PS (t] = 1) by about two orders of magnitude appears to be natural, because 
PAMS chains differ from those of PS by the presence of methyl groups attached to each tertiary 
backbone carbon [To] . 



6.2 Fitting of observations on poly(L— lactide) 

We proceed with the analysis of experimental data on thermal degradation of PLLA reported 
by Yu et al. [22]. PLLA samples with the initial mass-average molecular weight M w (0) = 108 
kg/mol and the polydispersity index D = 1.5 were heated to a required temperature T (in the 
range between 180 to 220 °C), preserved at this temperature in a reaction vessel under nitrogen 
atmosphere for given amounts of time (from 0.5 to 6.0 h), subsequently removed from the vessel, 
and their molecular weights were measured by GPC The ratio of the number- average molecular 
weights d n is plotted versus exposure time t in Figure 11. This figure shows that the time- 
temperature superposition principle is not valid for the degradation process (the experimental 
data at T = 200 and T = 220 °C re-plotted in the semi-logarithmic scale cannot be superposed 
by shifts along the time-axis with an acceptable level of accuracy). This conclusion may be 
explained by the fact that the activation energies E 1 and Er noticeably differ from each other. 
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To validate this result and to assess the difference between the activation energies for frag- 
mentation and annihilation, we approximate each curve depicted in Figure 11 separately. To 
reduce the number of adjustable parameters, we set a = (which means that the fragmentation 
rate is assumed to be independent of chains' length) and b = 2.5 (this value is taken because 
the observations at T = 200 and T = 220 °C depicted in Figure 11 resemble our results of 
numerical simulation with b = 2.5 reported in Figure 5). The fact that the rate of fragmen- 
tation in PLLA specimens with relatively small molecular weights (less than 300 kg/mol) is 
independent of chains' length is confirmed by the experimental data reported by Bywater and 
Black (see Figure 3 in [32 ). 

The rates 70 and T are determined by using an algorithm similar to that described in 
Section 5. We fix some intervals [0,7 max ] and [0, r max ], where the "best-fit" parameters 70 
and r are assumed to be located, and divide these intervals into J subintervals by the points 
7W = iA 7 and = jAT (i, j = 1, . . . , J - 1) with A7 = 7 max /J and Ar = r max /J. For any 
pair {7^,r^}, Eq. (35) with initial condition (33) is integrated numerically by the Runge- 
Kutta method with K = 100 and the time-step At = 0.1. The best-fit parameters 70 and To 
are determined from the condition of minimum of the cost function (46). Figure 11 reveals 
good agreement between the observations at all three temperatures T under consideration and 
the results of numerical simulation. The ratio rj determined by Eq. (34) ranges from 0.9 at 
T = 180 to 11.0 at T = 200 and to 15.0 at T = 220 °C. These values together with the results 
depicted in Figures 1 and 2 show that the effect of annihilation of end- and side-groups on 
the degradation process is quite substantial, and the influence of this process increases with 
temperature. 

The adjustable parameters 70 and r are plotted versus temperature T in Figure 12. The 
experimental data are approximated by the equations similar to Eq. (44), 

\nA = A Q -^, i nB = B -^ (47) 

with T ref = 453 K (180 °C) and 

^ = ~~ ref ' 

70 

B = pref ' 

Equations (47) follow from Eq. (42). Figure 12 shows that Eq. (47) provides an acceptable fit 
of the observations (it should be noted some scatter of the experimental data). However, the 
deviations of the data from Eq. (42) are noticeably less pronounced that the difference between 
curves 1 and 2 in Figure 12, which confirms our hypothesis that fragmentation and annihilation 
of chains may be governed by different processes at the micro-level. 

The activation energy found for the fragmentation process, = 91.9 kJ/mol, is in good 
agreement with the value of E 1 determined for PAMS. The fact that the activation energy 
for annihilation of end- and side-groups (Er = 223.9 kJ/mol) exceeds by about twice that for 
fragmentation of chains implies that the annihilation process strongly affects polymer degrada- 
tion at elevated temperatures, whereas its influence is of secondary impotence at relatively low 
temperatures. This conclusion provides a clue to explain our results of numerical simulation 
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(48) 



on PS specimens: the temperature T = 275 °C is relatively low for intensive breakage of end- 
and side-groups in polystyrene, which implies that the value 77 = 1.0 found by matching the 
experimental data is rather small. 



7 Thermo-gravimetric analysis 

It was revealed in the previous sections that the model can correctly describe experimental 
data on the evolution of molecular weights with exposure time in conventional depolymeriza- 
tion tests. It should be noted, however, that other kinetic models grounded on either the 
fragmentation concept or the aggregation-fragmentation theory can reproduce at least some 
of these observations with a reasonable level of accuracy. What cannot be described by the 
conventional models (all of which imply the conservation law for the number of segments in 
a network) is observations in thermo-gravimetric experiments, where a decrease in a sample's 
mass is measured as a function of exposure time at elevated temperatures. 

Our aim now is to demonstrate that (i) the model grounded on the fragmentation-annihilat- 
ion concept can qualitatively describe experimental data in thermo-gravimetric tests, and (ii) 
the kinetic equations with adjustable parameters determined by matching observations in GPC 
tests can quantitatively predict results of thermo-gravimetric analysis (TGA). 

Two kinds of TGA tests are conventionally performed. In experiments with a constant rate 
of heating q = dT/dt, mass m of a specimen is measured as a function of exposure time t and 
the ratio 

d(t) = -j^ 
m{0) 

of the current mass m(t) to the initial mass m(0) is plotted as a function of current temperature 
T(t). The graph d(T) demonstrates the following features: (i) d remains practically constant 
and equal unity up to some temperature T onset , at which the degradation process starts, (ii) d 
sharply decreases with temperature (within an interval of temperatures smaller than 100 K), 
and (iii) d practically vanishes at higher temperatures. The interval of pronounced changes in 
d is characterized by the depolymerization temperature T depol that describes the life-time of a 
polymer and equals the temperature at which the mass loss reaches 5 %. 

In GPA tests with a constant temperature T, a specimen is rapidly heated to the required 
temperature, which remains constant during the experimental procedure, and the sample's 
mass is measured as a function of exposure time t. The function d(t) monotonically decreases 
from d(0) = 1 and tends to zero with the growth of exposure time t. The decrease in mass at a 
given temperature T is characterized by the half-life time thi that is determined as the exposure 
time at which the sample loses half of its initial mass. 

To compare predictions of the model with the experimental data available in the literature, 
we perform numerical integration of Eqs. (33) and (35) with K = 100 and At = 1.0 -10~ 3 [these 
values guarantee that the evolution of the first moments M m (t) (m = 0,1,2) is independent 
of the integration algorithm] for poly(a;-methylstyrene) with a = 0.0, b = 1.5 and r\ = 50.0 
(these values are found by matching the experimental data presented in Figure 9). The effect 
of temperature T on the rate of fragmentation 70 is described by Eqs. (44) and (45) with 
7o ef = 9.0 • 10 -7 , A = 26.86 and A\ = 10344.0 (these values are determined by matching the 
observations depicted in Figure 10). To take into account the influence of the initial molecular 
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weight M w (0), we use Eq. (41) with M^ ef (0) = 9.0 kg/mol. The ratio of the current mass m(t) 
to the initial mass m(0) is determined by Eq. (31), where the moment Mi(i) is given by Eq. 
(12). 

Integration of Eq. (35) is carried out for M w (0) = 4.0 kg/mol and M w = 680.0 kg/mol 
(in the latter case, a correction factor of 1.5 is introduced in Eq. (41) in accord with the 
experimental data reported in Figure 3 of jH2j) for the first kind of thermo-gravimetric tests 
with q = 10.0 K/min and for M w = 680.0 kg/mol for the other kind of experiments at the 
temperature T = 560 K. The results of numerical analysis are presented in Figures 13 and 
14. These figures demonstrate that the shapes of the curves d(T) and d(t) are similar to those 
reported in numerous experimental studies. 

To show that our results of numerical simulation provide quantitative coincidence with 
observations, we depict (in Figure 13) the onset temperature T onset for degradation of PAMS 
with the low molecular weight M w (0) = 4.0 kg/mol [HE], and the depolymerization temperature 
Tdepoi for PAMS with the high molecular weight M w = 680.0 kg/mol [36J. For the same purpose, 
we plot the boundaries of the interval for the half-life time t^i of PAMS at T = 560 K reported in 
Striking similarity is observed between the experimental data and the model predictions. 

8 Conclusions 

A kinetic model is analyzed for thermal degradation of polymers. The governing relations are 
based on the fragmentation-annihilation concept. According to this approach, thermal degra- 
dation reflects two processes at the micro-level: (i) binary fragmentation of macromolecules and 
(ii) scission of end- and side-groups and their annihilation (diffusion and evaporation) from the 
network of chains. The evolution equations involve four material constants that are determined 
by fitting experimental data for polystyrene, poly(a-methylstyrene) and poly(L-lactide). 

Two explicit solutions of the integro-differential equations are found by using the method 
of generating function and the Charlesby method. These solutions (corresponding to special 
values of adjustable parameters) serve as a basis for numerical analysis. Numerical simulation 
of the kinetic equations is performed to study the effect of material constants on the evolution 
of mass-average and number-average molecular weights of polymers with exposure time. 

With reference to available experimental data, two superposition principles (time-temperat- 
ure and time-molecular weight) are formulated, and restrictions on material constants are found 
that guarantee their validity. These limitations allow the number of adjustable parameters to 
be reduced in order to make the approximation of observations more reliable. 

It is revealed that the fragmentation-annihilation concept not only can be applied to de- 
scribe experimental data on changes in the molecular weight, but can also be used to predict 
(qualitatively, and in some cases quantitatively) the mass decrease observed in TGA tests. 

The following conclusions may be drawn from our analysis of experimental data on PS, 
PAMS and PLLA: 

• The conventional assumption about the effect of chains' length (a > 0) on the fragmenta- 
tion process is excessive. It results in a noticeable complication of the theoretical analysis, 
but does not improve the quality of fitting. 

• The account of the influence of chains' length on the rate of annihilation is essential. The 
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parameter b weakly affects changes in the molecular weight with time at b e (0, 1), but 
its effect becomes important at higher values of the exponent b. 

• The kinetics of polymer degradation is noticeably affected by the annihilation process 
when two conditions are fulfilled: (i) the number of end- and side-groups in polymer chains 
is rather large, and (ii) the rate of annihilation exceeds the rate of fragmentation of chains 
by (at least) an order of magnitude. Our results of numerical simulation demonstrate that 
the latter condition is satisfied at relatively high temperatures. 



Appendix 

The aim of this section is to derive an explicit formula for the evolution of the second moment 
M 2 with time t, when the fragmentation process is described by the conventional equation 

dn f°° 

— (t,x) = -i n(t,x) + 270 n{t,y)dy (A-1) 

with an arbitrary initial condition (11). Equation (A-1) coincides with Eq. (11) where the 
annihilation process is disregarded (r = 0). The derivation is presented as an appendix, 
because we use a conventional approach (the method of generating functions) for the study of 
Eq. (A-1). To the best of our knowledge, the final result is, however, novel. 

We begin with the analysis of the moments M (t) and Mi(t). It follows from Eq. (15) that 
these functions are governed by the differential equation 

dMn , . „ , / s dMi . . . . . 

-5r<«> -ir^ = - < A - 2 > 

Integrating Eq. (A-2) and bearing in mind that 

M (0) = 1 

[this equality follows from Eqs. (1), (8) and (12)], we obtain 

M„(f) = 1 + 7oM 1 (0)t, Mx(t) = Mx(0). (A-3) 

We now multiply Eq. (A-1) by exp(— zx), where 2 is a new variable, and integrate the result 
over x from zero to infinity. Introducing the notation [the Laplace transform of the function 

n(t,x)\ 

n(t,z) = / n(t, x) exp(— zx)dx, (A-4) 



we find that 

— (t,z) — — 70 / xn(t, x) exp(—zx)dx + 270 / exp(— zx)dx \ n(t,y)dy. (A-5) 
at Jo Jo Jx 

According to Eq. (A-4), the first integral on the right-hand side of Eq. (A-5) reads 

/ xnit, x) exp(— zxjdx = — — — (t, z). 
Jo oz 
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To transform the other integral, we change the order of integration and use Eqs. (12) and 
(A-4), 



1 

z Jo 



roo roo roo ry 

/ exp(— zx)dx I n(t,y)dy = / n(t,y)dy / exp(— zx)dx 

JO Jx JO Jo 

1 - exp(-zy) n(t, y)dy = - M (t) - n(t, z) . 

L J % L J 

Substitution of these expressions into Eq. (A-5) results in the partial differential equation 

— (t, z) = lQ —(t, + [Mo(t) - n(t, z)\ . (A-6) 

Introducing the function 



$(t,z) =n(t,z) - M (t), 



(A-7) 



we present Eq. (A-6) in the form 



_( M ) + _ (t) = 70 _( M )-^ (t , 2 ). 



It follows from this equality and Eq. (A-2) that 

-(M)- 7o -(M) = -— 



$(t,z) -7oMi(0). 



We search a solution of Eq. (A-8) in the form $ = \J/(x!,x 2 ), where 



xi = t H , x 2 = z. 

7o 



(A-8) 



(A-9) 



In the new notation, Eq. (A-8) reads 

0* 



. (x 1: x 2 ) = M 1 (0) + — f(x u x 2 ) 
c*x 2 x 2 



(A-10) 



An advantage of Eq. (A-10) is that it includes xi as a parameter. We now introduce the 
function ^i(xi,x 2 ) by the formula 



\fr(xi,x 2 ) = x 2 *i(xi,x 2 ). 
Equations (A-10) and (A-ll) imply that this function satisfies the equation 

\Xi,X2) — • 



dx 



The general solution of Eq. (A-12) reads 



x 2 



(A-ll) 



(A-12) 



x 2 
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where F is an arbitrary smooth function of X\. Returning to the initial notation with the help 
of Eqs. (A-7), (A-9) and (A-ll), we obtain 



n(t, z) = M (t) - M 1 (0)z + Fit + —) 

7o 7 



It follows from Eqs. (11) and (A-4) that 

roo 

n(0, z) — no(z) : ho(z) — / exp(— zx)no(x)dx. 

Jo 

Setting t — in Eq. (A-13) and using Eq. (A-14), we find that 



F(z) 



1 



(loz) 2 

Equations (12) and (A-4) imply that 



n (7o^)-M (0) + M 1 (0)7o2 ; 



M 2 (t) 



d 2 h 



(t,0). 



Differentiation of Eq. (A-13) with respect to z results in 

d 2 n 



dz 2 



(t,0)=2F(t). 



Combining this equality with Eqs. (A-15) and (A-16), we arrive at the formula 

2 



M 2 (t) 



n (7ot)-M (0) + M 1 (0)7ot 



(A-13) 



(A-14) 



(A-15) 



(A-16) 



(A-17) 



Equations (A-3) and (A-17) provide explicit formulas for the evolution of three first moments 
M m driven by thermal degradation of a polymer with an arbitrary initial distribution of chain 
lengths. 
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List of figures 

Figure 1: The ratio of number- average molecular weights d n versus dimensionless time t. 
Curves: results of numerical simulation for a = 0.0, b = 1.0 and rj = 0.0, 1.0, 10.0, 100.0 
and 200.0, from top to bottom, respectively 

Figure 2: The ratio of mass- average molecular weights d w versus dimensionless time t. 
Curves: results of numerical simulation for a = 0.0, b = 1.0 and rj = 0.0, 1.0, 10.0, 
100.0 and 200.0, from top to bottom, respectively 

Figure 3: The ratio of number-average molecular weights d n versus dimensionless time t. 
Curves: results of numerical simulation for b = 1.0, rj = 10.0 and a = 0.3, 0.6, 0.9 and 
1.2, from top to bottom, respectively 

Figure 4: The ratio of mass-average molecular weights d w versus dimensionless time t. 
Curves: results of numerical simulation for b = 1.0, i] = 10.0 and a = 0.3, 0.6, 0.9 
and 1.2, from top to bottom, respectively 

Figure 5: The ratio of number-average molecular weights d n versus dimensionless time t. 
Curves: results of numerical simulation for a = 0.2, rj = 10.0 and b = 0.0, 1.0, 1.5, 2.0 
and 2.5, from top to bottom, respectively 

Figure 6: The ratio of mass-average molecular weights d w versus dimensionless time t. 
Curves: results of numerical simulation for a = 0.2, rj = 10.0 and b = 0.0, 1.0, 1.5, 
2.0 and 2.5, from top to bottom, respectively 

Figure 7: The ratio of mass-average molecular weights d w versus time t min. Symbols: 
treatment of observations on PS with the initial molecular weights M w (0) = 26 (triangles), 
110 (stars), 210 (asterisks), 330 (filled circles) and 930 kg/mol (unfilled circles). Solid line: 
results of numerical simulation 

Figure 8: The shift factor A versus the ratio of mass-average molecular weights. Circles: 
treatment of observations on PS. Solid line: approximation of the experimental data by 
Eq. (40) with A = 0.05 and A 1 = 0.83 

Figure 9: The ratio of number-average molecular weights d n versus time t. Symbols: treat- 
ment of observations on PAMS at the temperatures T = 110 (unfilled circles), T = 120 
(filled circles), T = 130 (asterisks), T = 140 (stars) and T = 150 °C (diamonds). Solid 
line: results of numerical simulation 

Figure 10: The shift factor A versus temperature T. Circles: treatment of observations on 
PAMS. Solid line: approximation of the experimental data by Eq. (44) with A = 26.86 
and A 1 = 1.03 • 10 4 

Figure 11: The ratio of number-average molecular weights d n versus time t. Symbols: ex- 
perimental data on PLLA. Solid lines: results of numerical simulation 
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Figure 12: The shift factors A (unfilled circles) and B (filled circles) versus temperature T. 
Circles: treatment of observations on PLLA. Solid lines: approximation of the experi- 
mental data by Eq. (47). Curve 1: A = 24.53, A 1 = 1.11 • 10 4 . Curve 2: B = 59.90, 
B 1 = 2.69 • 10 4 

Figure 13: The ratio d of a sample's mass to its initial mass versus temperature T. Solid 
lines: results of numerical simulation for PAMS under heating with the rate q = 10.0 
K/min. Curve 1: M w (0) = 4.0 kg/mol. Curve 2: M w (0) = 680.0 kg/mol. Vertical lines: 
the onset temperature T onset for thermal degradation of the low-molecular weight PAMS 
and the depolymerization temperature Td epo i for the high-molecular weight PAMS 

Figure 14: The ratio d of a sample's mass to its initial mass versus exposure time t at the 
temperature T = 560 K. Solid line: results of numerical simulation for PAMS. Vertical 
lines denote the interval of time where a specimen loses a half of its initial mass 
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